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ABSTRACT 


CLYDE  is  a computer  language  for  jrour  differ- 
ential equations.  It  provides  numerical  solutions  to 
an  important  class  of  second  order , elliptic  partial 
differential  equations  (Laplace  and  Poisson)  which  appear 
in  almost  every  branch  of  applied  mechanics:  governing  the 
solutions  to  design  problems  in  heat  conduction,  stress 
concentration,  and  potential  fields  (electric,  magnetic, 
electrostatic,  gravitation,  irrotational  fluid  flow,  etc..). 
There  are  three  versions  of  CLYDE.  This  document  describes 
the  capabilities  of  the  CDC  6000/TEKTRONIX  4014  storage 
tube  graphics  version  (CLYDE-TEK)  and  the  batch  version 
(CLYDE-B)  and  also  serves  as  a preliminary  user's  manual. 

An  earlier  version  (CLYDE-274),  written  for  the  CDC 
6500/1700/274  refresh  graphics  facility,  is  described  in 
MISD  Information  Report  73-1,  January  1973,  and  includes 
the  extension  of  the  solution  to  the  fourth  order  stress 
analysis  equation  for  flat  plates.  All  CLYDE  versions 
were  written  for  CDC  6000  host  computers  under  the  SCOPE 
operating  system  with  overlay  structures.  Two  applications 
covered  in  detail  in  this  document  are  steady  state  heat 
conduction  and  the  membrane  or  soap  film  analoqy  of  torsion 
of  bars  and  shafts. 
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BACKGROUND 


Most  armament  design  is  governed  by  the  classical 
ideas  and  equations  of  continuum  mechanics.  Throuqh  the 
descriptive  differential  equations  of  elasticity,  clas- 
sical mechanics,  electromagnetic  theory,  fluid  mechanics, 
etc.,  the  working  state  of  armament  items  may  be 
accurately  studied.  Physical  phenomena  in  continuous 
systems  - elastic  bodies,  fluids  - are  usually  described 
by  partial  differential  equations  with  their  associated 
boundary  or  initial  conditions.  Closed  form  solutions 
to  these  PDF's  are  rarely  available  in  the  design  room 
with  its  configurations  that  perversely  do  not  conform 
to  classical  text  book  illustrations.  Therefore,  in 
the  harsh  world  of  reality,  recourse  to  numerical 
solutions  is  an  absolute  necessity. 

One  widely  used  numerical  technique  is  finite 
differences.  CLYDE  solves  two  dimensional  boundary 
value  problems  with  generalized  contours,  using  finite 
difference  approximations.  A boundary  value  problem 
is  one  in  which  some  function (s)  of  the  problem 
variable  is  known,  but  only  at  the  boundary  of  the 
problem.  Steady  state  temperature  distribution  by 
means  of  heat  conduction  is  a good  example  of  this:  the 
temperatures  around  the  boundary  or  periphery  of  an 
arbitrary  shape  are  kept  constant  and  the  problem  is  to 
find  the  temperature  distribution  within  the  area  of  the 
shape . 


This  document  describes  two  versions  (the  latest 
interactive  graphics  and  the  basic  batch  versions)  of  a 
numerical  solution  to  an  important  class  of  boundary 
value  problems  involving  the  second  order  elliptic 
partial  differential  equations: 

= g(x.y) 


A_a£t  + B_9£f_ 
ax2  av2 


at 

3R 


9 (r,z) 
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To  illustrate,  consider  the  cross-section  of  an 
experimental  twin-lead  cable.  The  distribution  of  elec- 
trical potential  (E,  in  volts)  is  described  by  the  second 
order  partial  differential  equation  (called  a harmonic 
equation) : 

32E  . 0 
8XZ  8Y2 

The  problem  was  run  with  the  surface  of  one  lead  maintained 
at  +24  V.,  the  other  lead  at  -24V.,  while  the  outer  sheathing 
of  the  cable  was  kept  at  0 V.  The  resulting  computer  solution 
generated  a plot  showing  the  lines  of  constant  electrical 
potential.  The  finite  difference  grid  nodes  are  displayed  over 
only  half  the  cable  cross-section.  Since  the  problem  possessed 
symmetry,  it  was  only  necessary  to  work  with  the  smallest 
repeating  section  - in  this  case,  a semi-circle. 


As  another  example  of  CLYDE'S  application,  this  time  in 
cylindrical  coordinates,  consider  the  problem  of  stress 
concentration  in  a stepped  and  grooved  shaft  loaded  in 
torsion.  This  is  the  case  of  a solid  circular  cylinder  with 
varying  diameter  (the  collar  and  grooves).  The  governing 
compatibility  equation  in  terms  of  Saint- Venant *s  stress 
function  f(R,Z)  is: 


+ B 


( c at 

8R2  R SR 


g(r.z) 


where  g(r,z)*0  within  the  shaft  (r,z  domain)  and  some  con- 
stant value (s)  on  the  various  boundaries. 


At  the  far  ends  of  the  shaft,  away  from  the  discontin- 
uities of  contour,  the  stress  function  is  a function  of  the 
radius  only  - satisfying  the  equation: 

o)2f  5 cjf__Q 

oJr2  R dR 

The  abrupt  discontinuities  of  diameter  at  the  collar, 
however,  are  stress  concentrators  and  a judicious  solution 
of  the  governing  PDE  yields  visual  insight  of  the  stress 
pattern. 

Plots  of  the  different  values  of  the  stress  function 
represent  surfaces  of  revolution,  and  adjacent  plots  repre- 
sent the  inner  and  outer  surfaces  of  torque  transmitting 
tubes.  Plots  of  equally  spaced  values  of  ( Af*  constant) 
depict  equimomental  tubes  - that  is,  tubes  of  equal  torque  or 
moment  transmitting  capacity.  It  would  appear,  intuitively, 
that  the  thinnest  parts  of  these  tubes  represent  the  highest 
stressed  portions.  The  total  torque  transmitted  by  the  shaft 
is  the  sum  of  the  torque  transmitted  by  all  the  tubes. 

Another  application  is  the  charting  of  the  distortion 
of  the  earth's  electrostatic  field  by  terrain,  buildings  and 
even  earthlings.  Both  examples  are  covered  in  more  detail 
later. 
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PROCEDURE 

The  resultinq  computer  program (a)  called  CLYDE 
(Computer  Lanquaqe  for  Your  Differential  Equations) 
employs  those  mathematical  models  already  in  the  finite 
difference  literature.  What  is  different  is  the  unique 
blending  of  these  models  with  ARRADCOM  developed  pattern 
recognition  algorithms  in  a man-oriented  input  and  output 
environment  of  a computer  time-sharing  facility. 

The  picture  or  contour  (mathematically  speaking,  the 
enclosed  domain)  of  the  problem  is  inputted  to  the  program, 
generally  on  punched  cards,  as  a series  of  contiguous 
straight  lines  and  circular  arc  segments.  The  computer 
program  overlays  the  picture  (or  domain)  of  the  problem 
with  a rectangular  network  of  vertical  and  horizontal  grid 
lines  - the  finite  difference  grid.  To  conserve  computer 
core  and  buffer  storage  space,  the  graphics  program  does 
not  always  display  the  complete  grid  lines  on  the  screen. 
Usually,  only  the  intersections  of  these  grid  lines  are 
displayed,  shown  as  little  crosses  or  plus  signs.  Aqain, 
to  conserve  computer  storage,  a mirror  line  algorithm  was 
developed  to  enable  users  to  take  advantage  of  any 
symmetry  of  the  problem,  and  analyze  only  a "repeating 
section"  of  the  problem  domain.  The  intersections  of 
the  horizontal  and  vertical  grid  lines  with  the  boundary 
are  called  boundary  nodes  and  are  shown  in  the  graphics 
version  as  little  circles.  The  intersections  of  the 
horizontal  and  vertical  grid  lines  (little  crosses) 
within  the  closed  boundary  of  the  problem  are  known  as 
inner  domain  nodes,  and  it  is  here  that  the  solution  is 
desired. 

At  each  inner  domain  node,  a finite  difference  ap- 
proximation to  the  governing  partial  differential  equa- 
tion (PDE)  is  generated  by  CLYDE.  The  resulting  set  of 
linear  algebraic  equations  are  solved  simultaneously  by 
the  program  for  the  unknown  problem  variable  (temper- 
ature, voltage,  stress  function,  etc..)  at  each  node  in 
the  overlaying  finite  difference  qrid.  The  graphics 
version  also  generates,  and  displays  on  the  screen, 
iso-value  contour  maps  for  any  desired  values  of  the 
variable  just  solved  for.  This  way,  a more  meaningful 
picture  of  the  solution  is  made  available  to  the  engi- 
neer; in  the  form  of  temperature  distributions,  constant 
voltage  lines,  stress  concentration  graphs,  or  even 
contour  lines  of  different  values  of  deformation  and 
bending  moment  in  structural  problems. 
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The  user  may  also  specify  a finer  grid  spacing  to  in- 
crease resolution  in  critical  regions  of  the  problem,  modify 
the  scale  of  the  display,  change  the  boundary  of  the  problem 
(or  redraw  it  completely) , and  change  boundary  conditions  and 

coefficients  all  at  the  face  of  the  screen.  It  is  also 

possible  to  request  CLYDE  to  pass  a plane  through  the  two 
dimensional  picture  displayed  on  the  screen.  This  plane 
is  perpendicular  to  the  screen  and  shows  as  a straight  line. 
CLYDE  will  generate  a new  display  showing  a cross  section 
(or  elevation)  view  from  the  edge  or  side.  In  this  manner 
the  variation  or  plot  of  the  solved  variable  along  that  line 
is  displayed  on  the  screen.  If  the  problem  geometry  is 
symmetrical,  the  designer  does  not  have  to  display  and  work 
with  the  entire  picture  of  the  problem.  If  he  desires,  he 
need  only  to  display  the  "repeating  section".  (This  is  done 
in  the  illustrative  solutions) . In  essence,  the  graphics  user 
may  examine  the  problem  solution  at  will,  and  redesign  the 
problem  at  the  screen  (problem  contour,  boundary  conditions, 
equation  coefficients,  etc.)  and  resolve  the  "new  design- 
problem. 

I 

\ 
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ENGINEERS ' QUICK  LOOK  AT  FINITE  DIFFERENCES 


The  finite  difference  approximations  to  the  partial 
differential  operators,  substituted  into  the  governing 
partial  differential  equation  (PDE) , yield  an  algebraic 
approximation  equation.  One  such  algebraic  equation  is 
generated  by  the  computer  program  at  each  INNER  DOMAIN 
NODE.  For  example,  substituting 


an., 

ay* 


H-2W 


into  the  equation 


+ B 


ait 

8Y2 


g(x.y) 


D 


and  letting  hx  = hy  = h (for  a square  grid)  yields  the 
following  algebraic  expression  called  the  harmonic 
operator  at  a typical  node  (labeled  node  0) : 


A (f j+f 3) +B (f 2+f 4 ) - (A+B) 2f 0=h2D 


This  finite  difference  equation  at  node  zero  involves 
the  unknown  variable  at  node  zero  (fo)  plus  the  unknown 
value  of  the  variable  at  the  four  surrounding  nodes 
(fl, 3*^4) > plus  the  grid  spacing  (h) . The  five  nodes 
involved  form  a four  arm  star  with  node  zero  at  the  center. 
This  algebraic  or  difference  equation  could  be  conveniently 
visualized  as  a four  arm  computation  stencil  made  up  of 
five  "balloons"  connected  in  this  four  arm  star  pattern 
and  overlayed  on  the  grid  nodes.  The  value  within  each 
balloon  is  the  coefficient  by  which  the  variable  (f)  at 
that  node  is  multiplied  to  make  up  the  algebraic  approxima- 
tion equation. 


# (Complete  mathematical  why's  and  how's  in  Appendix  b) 
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USING  CENTRAL  DIFFERENCES, THE  FINITE 
DIFFERENCE  APPROXIMATIONS  TO  THE  PARTIAL 
DIFFERENTIAL  OPERATORS  , OF  THE  FUNCTION 
f,  AT  REPRESENTATIVE 


NODE  0 ARE  : 

’"ty-  * ZHT 

■glL.  _L_(f,-2f„*f5) 

0 x*  h; 

£.lf . s — 1— . ( f - ? f ♦ f i 

ay*  h*  v * 0 v 


PROGRAM  DATA 


INPUT  DATA  CARD  PREPARATION 

The  complete  set  of  input  data  and  production  options 
must  be  inputted  to  the  CLYDE-BATCH  program  on  a deck  of 
(up  to  nine)  different  types  of  punched  cards.  Only 
four  of  these  nine  types  are  required  to  initialize  the 
interactive  graphics  CLYDE-TEK  version,  since  problem 
customizing  is  performed  and  output  options  selected 
at  the  graphics  screen.  The  subset  of  four  card  types 
required  for  CLYDE-TEK  initialization  are  noted  with  a 
"(TEK)"  typed  to  the  right  of  the  card  type. 


The  nine  types  of  input  data  cards  are: 

Type  1 IDENTIFICATION  CARD  (TEK) 

Type  2 OPTION  CARD 

Type  3 PROBLEM  DEFINITION  CARD  (TFK) 

CONTOUR  LINE  SEGMENT  IDENTIFICATION  CARD,  sets  of  2 


Type  4 LINE  SEGMENT 

Type  5 BOUNDARY  CONDITIONS 
Type  6 FINER  GRID  CARDS 

PLOT  INSTRUCTION  CARDS,  sets  of  3 
Type  7 PLOT  CONTROL  CARD 

Type  8 CONTOUR  PLOT  CARD 

Type  9 CROSS  SECTION  PLOT  CARD 


(TEK) 

(TEK) 
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(Card  Type  1)  IDENTIFICATION  CARD  (TEK) 
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GENERAL  NOTES  ON  INPUT  DATA 


CONTOUR  LINE  SEGMENT  CARDS  - The  boundaries  of  the 
problem  (the  outer  and  all  inner  contours)  are  inputted  to 
CLYDE  as  polystrinqs  of  contiguous  straight  lines  and  cir- 
cular arc  seqments.  The  line  segment  sequence  for  any  con- 
tour or  boundary  may  beqin  with  any  line  seqment  of  that 
contour  and  may  proceed  in  either  a clockwise  or  counter- 
clockwise direction  around  that  boundary.  Once  begun,  how- 
ever, the  segments  must  be  input  in  order,  in  the  direction 
started,  until  that  boundary  is  completed.  Each  LINE  SEG- 
MENT card  must  be  followed  by  a BOUNDARY  CONDITIONS  card. 
There  are  NCPT  pairs  of  these  (LINE  SEGMENT  and  BOUNDARY 
CONDITIONS)  cards. 

Only  the  circular  arcs  (CA)  and  equation  arcs  (EA)  re- 
quire origin  coordinate  (XO,yO)  and  direction  of  development 
(DIR)  input  data.  These  fields  are  ignored  by  the  prooram 
when  the  line  segments  are  straight  lines  (SL) , mirror 
lines  (ML) , or  equation  lines  (EL) . 

A mirror  line  (ML)  is  a line  of  symmetry,  a veritable 
line  of  reflection.  Its  primary  function  is  to  permit  the 
analysis  of  a "repeating  section"  of  the  problem.  The  same 
number  of  qrid  lines  overlaying  this  smaller,  but  repre- 
sentative, portion  as  is  normally  used  over  the  entire 
problem  area  produces  a finer  finite  difference  net.  This 
yields  higher  resolution  and  a more  accurate  numerical 
approximation  while  requiring  no  more  storage  or  running 
time  - a delightful  freebie.  Unfortunately,  the  present 
state  of  CLYDE'S  development  allows  for  horizontal  and 
vertical  mirror  lines  only,  so  that  a quadrant  is  the  small- 
est symmetrical  portion  that  can  be  presently  handled.  The 
user  is  cautioned  to  position  his  problem  so  that  any 
vertical  mirror  line  is  the  leftmost  line  in  the  data  set 
and  any  horizontal  mirror  line  the  lowest  horizontal  line 
in  the  data  set.  Incidentally,  a mirror  line  also  implies 
(and  may  be  so  used)  a boundary  condition  of  the  first 
derivative  of  the  problem  variable,  normal  to  the  boundary, 
being  zero.  a^n-o 

The  equation  line  (EL)  and  equation  arc  (EA)  merely 
indicate  a linear  variation  in  the  boundary  value  of  the 
problem  variable  from  one  end  of  the  line  (or  arc)  segment 
to  the  other,  along  the  length  of  the  line  (EL)  or  the  curve 
of  the  arc  (EA) . 
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Tk*  total  number  of  contour  segments  in  any  problem 

may  not  exceed  125.  The  total  number  of  inner  and  outer 
contour  nodes  (intersections  of  contour  segments  and  grid 
lines)  may  not  exceed  500.  The  total  number  of  inner 
domain  nodes  (intersections  of  horizontal  and  vertical 
grid  lines  within  the  problem  area)  may  not  exceed  1000. 
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SOLUTION  MATRIX  fc  BAND  WIDTH 


CLYDE  generates  a finite  difference  equation  approxi- 
mation to  the  PDE  to  be  solved  at  each  inner  node  within 
the  problem  area  or  domain.  For  a relatively  straight 
forward  ten  unit  square  with  a unit  grid  this  would 
produce  81  equations  (one  for  each  of  the  9x9  or  81  nodes) 
with  81  unknowns,  to  be  solved  simultaneously.  Consider- 
ing a large  complex  problem,  with  a realistically  (useful) 
fine  grid,  the  number  of  equations  and  unknowns  becomes 
horrendously  formidable.  Realizing,  however,  that  there 
are,  at  most,  five  non-zero  coefficients  in  any  of  the 
equations  it  was  decided  to  keep  the  non-zero  terms 
clustered  "tightly”  about  the  main  diagonal  and  use  a linear 
equation  solution  routine.  This  routine  provides  a "neat" 
computer  procedure  for  solving  our  anticipated  N X N system 
of  linear  simultaneous  equations  whose  coefficient  matrix 
is  of  the  band  form  (i.e.,  it  has  non-zero  elements  only 
about  the  main  diagonal  and  zeroes  elsewhere) . Only  the 
band  elements  need  be  stored,  permitting  the  solution 
of  large  systems  of  linear  simultaneous  equations  in 
relatively  few  storage  locations.  Gaussian  elimination  is 
used,  modified  to  take  advantage  of  the  reduced  matrix. 

The  routine  also  uses  partial  pivoting  to  reduce  roundoff 
error. 

The  maximum  band  width  programmed  into  the  present 
CLYDE  is  90  and  exceeding  this  will  produce  a program  stop 
and  a message  as  shown  on  the  accompanying  illustration. 

A box  is  displayed  about  the  node  whose  total  bandwidth 
exceeds  90  or  about  the  node  whose  "left  hand"  or  "right 
hand  bandwidth"  exceeds  45.  The  band  width  at  any  node 
is  the  total  number  of  nodes  from  its  neighbor  immediately 
to  its  left  to  its  immediate  right  neighbor.  The  node 
count  is  up  the  vertical  grid  lines  from  the  left 
neighbor,  through  the  node  in  question,  and  up  to  the 
right  neighbor.  Boundary  nodes  are  not  neighbors,  nor 
are  they  counted.  The  CLYDE  user  should  be  carefully 
selective  in  his  choice  of  initial  grid  spacing, 
problem  orientation,  shape  and  size  of  the  finer  mesh 
if  any,  and  repeating  section  options  to  avoid  exceeding 
the  band  width  limit. 
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STEP-BY-STEP  ILLUSTRATIVE  EXAMPLES 
STEADY-STATE  TEMPERATURE  DISTRIBUTION  (X,Y  coordinates) 

Heat  conduction,  at  steady  state  conditions,  within 
a manifold  with  irregular  channels  is  examined  on  the 
following  pages.  The  square  cross-section  of  manifold  with 
its  circular  and  slotted  perforations,  possesses  quadrant 
symmetry.  It  is  this  quarter  of  cross-section  area  that 
is  input  to  the  CLYDE  program  in  cartesian  coordinates. 

The  outer  perimeter  of  the  square  is  maintained  at  a 
constant  temperature  T^  of  0 P.  The  contours  of  the  two 
inner  circular  holes  are  maintained  at  a constant  tempera- 
ture Tj  of  50*  F,  while  the  boundary  of  the  double-hole- 
slot  is  kept  at  yet  a third  constant  temperature  T.  of 
100  F.  The  problem  is  to  determine  the  temperature  dis- 
tribution throughout  the  cross-section  of  the  manifold. 
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STEADY-STATE  TEMPERATURE  DISTRIBUTION  (R,Z  coordinates) 

It  is  necessary  to  know  the  steady-state  temperature 
distribution  within  the  walls  of  an  experimental  convergent 
nozzle  of  a sled  rocket.  The  ambient  temperature  at  the 
outer  skin  of  the  nozzle  is  a constant  72*F.  The  inner  * 
walls  of  the  nozzle  exhibit  a uniform  temperature  of  1000  F 
at  the  combustion  chamber  and  linearly  decreasing  in  the 
converging  nozzle  portion  to  300*F  at  the  throat.  The  wall 
of  the  cylindrical  combustion  chamber  has  temperatures  vary- 
ing radially  from  1000*F  at  the  inside  to  72*F  at  the  outside. 
The  temperatures  within  the  wall  at  the  throat  vary  from 
300*F  on  the  inside  to  72*F  at  the  outside  skin.  The  prob- 
lem, similar  to  that  of  the  manifold,  is  to  calculate  and 
plot  the  thermal  profiles  within  the  walls  of  the  chamber- 
nozzle.  In  cylindrical  coordinates,  the  Z replaces  the 
X axis  of  the  previous  problem,  while  R replaces  the  Y. 
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PUNCH  CARD  INPUT  FOR  CLYDE-TEK 
RUN  OF  NOZZLE  THERMAL  PROBLEM. 
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ax2  9Y2 


Consider  the  torsion  of  a 1 inch  diameter  shaft  cut 
by  two  diametrically  opposed  identical  key  ways.  The 
symmetry  of  this  prismatic  bar  permits  the  investigation 
of  only  one  quadrant.  The  cross-section  was  enlarged  by 
a factor  of  10  (no  effect  on  results)  and  all  dimensions 
(X,Y  values)  are  in  terms  of  these  magnified  units.  The 
governing  equation  is  shown  above  and  the  boundary  con- 
dition is  O* 
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CARD  INPUT  FOR  CLYDE-TEK 
QUADRANT  OF  2 -KEYWAY 


(Z)  AND  (X)  KEYS  PRODUCE  CROSS- 
SECTION  OR  SIDE  VIEWS  OF  STRESS 
FUNCTION  VARIATION. 
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PUNCH  CARD  INPUT  FOR  TORSION 
OF  TRIANGULAR  BAR. 
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VARIATION  AT  MIRROR  LINE 
DISPLAYING  PULL  VIEW. 
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CLYDE- TEK.  TORSION  (EQUILATERAL  TRIANOLE  CROSS-SECTION)  OAR.  l/tB/77 
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*(See  example,  next  page) 


AUTO  DELETE:  Selecting  this  command  will  automatically  delete 

all  extraneous  nodes,  that  is,  all  domain  nodes 
outside  the  problems'  boundaries  as  well  as  all 
nodes  within  holes  or  inner  contours  of  the  prob- 
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1 produce  five  contour  maps  for  the  values  of 
, 150,  700,  250,  and  300. 
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RESTART  (NO  NODES) : Same  as  above,  but  picture  is  redrawn  with  NO 

NODES  displayed. 
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THE  TORQUE- ING  CHAPTER 


The  elastic  stress  analysis  of  uniformly  circular 
shafts  in  torsion  is  a familiar  and  straightforward 
concept  to  design  engineers.  As  the  bar  is  twisted, 
plane  sections  remain  plane,  radii  remain  straight, 
and  each  section  rotates  about  the  longitudinal  axis. 
The  shear  stress  at  any  point  is  proportional  to  the 
distance  from  the  center  and  the  stress  vector  lies 
in  the  plane  of  the  circular  section  and  is  perpen- 
dicular to  the  radius  to  the  point,  with  the  maximum 
stress  tangent  to  the  outer  face  of  the  bar  (another 
shearing  stress  of  equal  magnitude  acts  at  the  same 
point  in  the  longitudinal  direction) . The  torsional 
stiffness  is  a function  of  material  property,  angle 
of  twist,  and  the  polar  moment  of  inertia  of  the 
circular  cross-section.  These  relationships  are 
expressed  as: 


& - T/J*  G,  or  T 


and  S - T*r/J,  or  S 

S 9 


G*©*r 


where  T « twisting  moment  or  transmitted  torque;  G = 
Modulus  of  Rigidity  of  the  shaft  material;  ©=■  angle  of 
twist  per  unit  length  of  the  shaft;  J ■ polar  moment  of 
inertia  of  the  (circular)  cross-section;  Sg  * shear 
stress;  and  r * radius  to  any  point. 

However,  if  the  cross-section  of  the  bar  deviates 
even  slightly  from  a circle,  the  situation  changes 
radically  and  far  more  complex  design  equations  are 
required.  Now,  sections  of  the  bar  do  not  remain 
plane,  but  warp  into  surfaces,  and  radial  lines  through 
the  center  do  not  remain  straight.  The  distribution  of 
shear  stress  on  the  section  is  no  longer  linear  and  the 
direction  of  shear  stress  is  not  normal  to  a radius. 

The  governing  equation  of  continuity  (or  com- 
patibility) from  Saint-Venant ' s theory  is 


82^  . 82^ 
8X2  8Y2 


2G© 
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where (f)m  Saint-Venant ' s torsion  stress  function.  The 
problem  then  is  to  find  a $ function  which  satisfies 
tf»is  equation  and  also  the  boundary  conditions  that 
<P  ■ a constant  along  the  boundary.  This  <p  function 
has  the  nature  of  a potential  function,  such  as 
voltage,  hydrodynamic  velocity,  or  gravitational 
height.  Its  absolute  value  is,  therefore,  not  important; 
only  relative  values  or  differences  are  meaningful. 

The  solutions  to  this  equation  required  com- 
plicated mathematics  and  even  simple,  but  commonplace, 
practical  cross-sections  could  not  be  reduced  to 
mathematical  formulae  - and  numerical  approximations  or 
intuitive  methods  had  to  be  used. 

One  of  the  most  effective  numerical  methods  to 
solve  for  Saint-Venant' s torsion  stress  function  is 
that  of  finite  differences.  The  CLYDE  computer 
program  was  applied  to  a number  of  keywayed  shafts  to 
produce  the  dimensionless  design  charts  on  the 
following  pages.  Most  of  the  charts  required  approxi- 
mately 45  computer  runs  for  plot  data  generation  - 
i but  once  completed,  the  design  charts  for  that 

' cross-section  are  good  for  virtually  all  combinations 

. of  dimensions,  material,  and  shaft  twist. 


The  three  dimensional  plot  of  (k  over  the  cross- 
section  is  a surface  and,  with  (f>  set  to  zero  (a 
perfectly  valid  constant)  along  the  periphery,  the 
surface  is  a domb  or  (t  membrane.*  It  has  been  proven 
that  the  transmitted  torque  (T)  is  proportional  to 
twice  the  volume  under  the  membrane  and  the  stress 
(S  ) is  proportional  to  the  slope  of  the  membrane  in 
the  direction  perpendicular  to  the  measured  slope.  For 
bars  with  solid  cross-sections,  the  maximum  stress 
(neglecting  the  stress  concentration  of  sharp  re-entrant 
corners  which  are  relieved  with  generous  fillets)  is  at 
the  point  on  the  periphery  nearest  the  center. 


*The  best  intuitive  method,  incidentally,  came 
from  Prandtlt  the  membrane  analogy.  He  showed  that 
the  compatibility  equation  for  a twisted  bar  was  the 
"same”  as  the  equation  for  a membrane  stretched  over 
a hole  in  a flat  plate  and  then  inflated.  This 
concept  provides  a simple  way  to  visualize  the 
torsional  stress  characteristics  of  shafts  of  any 
cross-section  relative  to  those  of  circular  shafts 
for  which  an  (exact)  analytical  solution  is  readily 
obtainable. 
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TORSIONAL  PROPERTIES 
OF 

SOLID,  NON-CIRCULAR  SHAFTS 


T - TRANSMITTED  TORQUE  (IN  - LBS) 

0 - ANGLE  OF  TWIST  PER  UNIT  LENGTH  (RADS/IN) 

Q - MODULUS  OF  RIGIDITY  OR  MODULUS  OF 
ELASTICITY  IN  SHEAR  (LBS/IN  2) 

R - OUTER  RADIUS  OF  CROSS-SECTION  (IN) 

V,  4^  - VARIABLES  FROM  CHARTS  (OR  TABLES) 

RELATED  TO  VOLUME  UNDER  “SOAP  FILM 
MEMBRANE"  AND  SLOPE  OF  “MEMBRANE" 

8S  - SHEAR  STRESS  (LBS/IN2) 
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2.1 


SOME  COMMENTS  ON  THE  ACCURACY  OF  THE  SOLUTION 


A comparison  may  be  easily  made  between  the  results 
produced  by  a CLYDE  analysis  of  the  torsion  of  a solid 
circular  shaft  and  the  (exact)  classical  text  book  solution. 
One  quadrant  of  a unit  raduis  solid  circular  shaft  was  run 
with  two  finite  difference  grid  spacings  and  the  results 


of  the  following  sets 

of  equations  are 

compared: 

EQUATION  COMPARISON 

CLYDE 

EXACT 

TORQUE 

2G  & (V)  R4  - 

m +-  G&J 

2 (V)  R4 

- — J 

2 (V)  R4 

•-  dr/2)R4 

2V 

— ► (fl/2) 

SHEAR  STRESS  (MAX) 


i 

t 


1 

CLYDE 

EXACT 

% DEVIATION 

TORQUE 

(h-0.125) 

1.5546 

1.5708 

1.03 

(h-0.0625) 

1.5669 

1.5708 

0.25 

SHEAR  STRESS 

(h-0.125) 

0.9379 

1.0 

6.21 

(h-0.0625) 

0.9688 

1.0 

3.125 

AREA* 

(h-0.125) 

3.13316 

3.14159 

0.268 

(h-0.0625) 

3.13984 

3.14159 

0.056 

♦(Used  for  internal  program  checking) 
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ONE  QUADRANT  OF  THE  UNIT  RADIUS 
CIRCULAR  SHAFT  WITH  FINITE 
DIFFERENCE  GRID. 
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FUNCTION  VALUES  ADDED. 
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APPENDICIES 

A.  OTHER  APPLICATIONS 

B.  MATHEMATICAL  MODEL 

C.  THINGS  TO  COME 

D.  AUTHORS'  CONTROL  CARDS 

E.  REFERENCES 

F.  UPDATE  SUBSCRIPTION  SERVICE  CARD 
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OTHER  APPLICATIONS 


MONITORING  EARTH'S  ELECTROSTATIC  FIELD 


An  evolutionary  version  of  the  CLYDE  program  has 
been  used  to  investigate  the  feasibility  of  portable 
anti-intrusion  devices  and  autopilot  sensors.  Both 
are  based  upon  the  same  principle:  the  atmosphere 
being  a giant  capacitor.  The  earth  is  negatively  charged 
and  the  air  above  positively  charged  with  voltage  layers 
ranging  from  zero  at  the  surface  ground  to  about  350,000 
volts  at  the  top  of  the  atmospheric  layer.  The  voltage 
gradient  is  greatest  at  sea  level  (50  volts/foot)  and 
decreases  with  altitude.  The  current  is  too  small  to 
sense,  but  the  voltage  layers,  nevertheless,  do  exist, 
and  can  be  measured. 

Anti- Intrus ion  Detectors : 

A man  standing  in  an  open  field  distorts  the  earth's 
electrostatic  field.  Houses  and  boulders  also  distort  the 
field,  but  they  are  stationary  while  an  intruder  must 
normally  move  (to  intrude  effectively).  The  anti-intrusion 
device  is  designed  to  monitor  the  rate  of  field  distortion 
caused  by  the  moving  intruder. 

Autopilot  Sensor: 

Roll  and  pitch  sensors  mounted  on  wing  tips  and 

inear  nose  and  tail  of  the  aircraft  read  different  volt- 
age signals  when  flight  deviates  from  the  horizontal  plane. 
Pairs  of  these  signals  (wing  tips,  nose  and  tail)  fed  into 
differential  voltmeters  drive  the  roll  and  pitch  servos 
to  correct  the  flight.  The  system  works  well,  so  far,  in 
nonmetallic  models  in  good  weather. 
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CONTOUR  LINES  OP  CONSTANT 
(ELECTROSTATIC)  VOLTAGE.  ARMS 
AT  90°. 


PUNCHED  CARD  CLYDE-TEK  INPUT 
FOR  STEPPED  AND  GROOVED  CIRCULAR 
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QUADRANT  SYMMETRY  SIMPLIFIES 
DATA  PREPARATION. 
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FINAL  DESIGN  - SMOOTH  CURVES 
ELIMINATE  STRESS  CONCENTRATORS 


MATHEMATICAL  MODEL 


As  the  term  implies,  boundary  value  problems  are 
those  for  which  conditions  are  known  at  the  boundaries. 
These  conditions  may  be  the  value  of  the  problem 
variable  itself  (temperature,  for  example) , the  normal 
gradient  (or  variable  slope),  or  higher  derivatives  of 
the  problem  variable.  For  some  problems,  mixed  boundary 
conditions  may  have  to  be  specified:  different 
conditions  at  different  parts  of  the  boundary.  CLYDE 
solves  those  problems  for  which  the  problem  variable, 
itself,  is  known  at  the  boundary. 


Given  sets  of  equally  spaced  arguments  and 
corresponding  tables  of  function  values  the  finite 
difference  analyst  may  employ  forward,  central,  and 
backward  difference  operators.  CLYDE  is  based  upon 
the  central  difference  operators  to  approximate  each 
differential  operator  in  the  equation. 

The  problem  domain  is  overlaid  with  an  appropriately 
selected  grid.  There  are  many  shapes  (and  sizes)  of 
overlaying  cartesian  and  polar  coordinate  grids: 

rectangular. . 
square. . 

equilateral  - triangular. . 
equilangular  - hexagonal . . 
oblique. . 


CLYDE  uses  a constant  size  (throughout  the  area  of 
the  problem)  square  grid  for  which  the  percentage  errors 
are  of  the  grid  size  squared  (h  *)  .'  This  grid  or  net 
consists  of  parallel  vertical  lines  (spaced  h units 
apart)  and  parallel  horizontal  lines  (h  units  apart)  which 
blanket  the  problem  area  from  left-to-right  and  bottom- 
to-top. 


ft 
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The  intersection  of  the  grid  lines  with  the 
boundaries  of  the  domain  are  called  boundary  nodes. 

The  intersections  of  the  grid  lines  with  each  other 
within  the  problem  domain  are  called  inner  domain  nodes. 
It  is  at  these  inner  domain  nodes  that  the  finite 
difference  approximations  are  applied.  The  approximation 
of  the  partial  differential  equation  with  the  proper 
finite  difference  operators  replaces  the  PDE  with  a 
set  of  subsidiary  linear  algebraic  equations  - one  at 
each  inner  domain  node.  For  the  practical  application 
of  the  method,  it  must  be  capable  of  solving  problems 
whose  boundaries  may  be  curved.  In  such  cases, 
boundary  nodes  are  not  all  exactly  h units  away  from 
an  inner  node  as  is  the  case  between  adjacent  inner 
nodes.  The  finite  difference  approximation  (of  the 
harmonic  operator)  at  each  inner  node  involves  not 
only  the  variable  value  at  that  node  and  at  the  four 
surrounding  nodes  (above,  below,  left  and  right)  but 
also  the  distance  between  these  four  surrounding  nodes 
and  the  inner  node  - and  at  the  boundaries  these  dis- 
tances, quite  likely,  vary  unpredictably . Compensation 
for  the  variation  of  these  distances  must  be  included 
in  the  finite  difference  solution.  CLYDE  represents 
the  problem  variable  by  a second  degree  polynomial  in 
two  variables,  and  employs  a generalized  irregular  star 
in  all  directions  for  each  inner  node.  In  practice, 
one  should  not  select  such  a coarse  grid  that  more  than 
(or  even)  two  arms  of  the  star  are  irregular  (or  less 
than  h units  in  length) . The  generalized  star  permits 
(and  automatically  compensates  for)  a variation  in  length 
of  any  of  the  four  arms  radiating  from  a node.  For  no 
variation  in  any  arm,  the  algorithm  reduces  exactly  to 
the  standard  harmonic  "computation  stencil". 


CONSIDER  THE  GENERAL  EXPRESSION: 

V'f.A^+B^-t-ff.0  EQd) 

IN  THE  17. £,X  COORDINATE  SYSTEM, 

WHERE  A,  B,C,D  ARE  ARBITRARY  CONSTANTS. 

WHEN  C«0,  V2f  REDUCES  TO  A TWO  COORDINATE 
SYSTEM,  SAY  IN  X AND  Y: 

V2f«A^2+  B^2«D  . . EQ  (2) 

USING  CENTRAL  DIFFERENCES, THE  FINITE 
DIFFERENCE  APPROXIMATIONS  TO  THE  PARTIAL 
DIFFERENTIAL  OPERATORS,  OF  THE  FUNCTION 
f , AT  REPRESENTATIVE 
NODE  0 ARE  : 

$£*•557  *2KT 

-|X  = 4_(f,-2f,+  fs) 
d^x*  h* 
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I. HARMONIC  OPERATOR  FOR  SQUARE  GRID 


FOR  A SQUARE  GRID,hx«  hy  * h , AND  THE  HARMONIC 
OPERATOR  V2f  BECOMES  *. 
h*V2f0-  [ A ♦ B (ft  + f4)-  ( A+B)  2f0]«  h*  D 


. .EQ  (3) 


SEE  FIS . I 


THE  NUMERICAL  TREATMENT  OF  AN  IRREGULAR 
STAR  (h,*  hKPh,*h4)  REPRESENTS  THE  FUNCTION  f, 
NEAR  THE  REPRESENTATIVE  NODE  0,  BY  A SECOND 
DEGREE  POLYNOMIAL  IN  X AND  Y • 
f ( X ,Y ) * f0+  a,X*  a*Y  + as  X2+  a4YVi»XY 
EVALUATING  THIS  POLYNOMIAL  AT  THE  NEIGHBORING 
NODES  (1, 2, 3, 4)  PRODUCE  THE  FOLLOWING  SET  OF 
EQUATIONS: 
f,  8 f0  ♦ 3i  h,  + ash,f 
f*  * to  + o* 

fs  * fo  “ 0|  hj  ♦ O3 
f4  * to  - 0*  h4  ♦ a4  h4 


WHICH  ARE  THEN  SOLVED  FOR  05  AND  04, WHICH 
ARE  NECESSARY  TO  SATISFY  THE  HARMONIC 
OPERATOR  V*f,  SINCE  • 

•a,+  2a3X  + q5  Y,~  ^2«  2as 

* <»*  *2a4Y  + a8  X,  = 2a4 
O Y O' 

AND 

A ( 2 a3)  + B (2a4) 

PERFORMING  THE  NECESSARY  ALGEBRAIC  OPER- 
ATIONS, SUBSTITUTING  RESULTS, COLLECTING 
TERMS  ,AND  USING  THE  FOLLOWING  RATIOS 


IRREGULAR  STAR  AT  NODE  0 


NEIGHBORING  NODES  ( I,  2,3,4, ) 
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THE  HARMONIC  OPERATOR  BECOMES : 


h2V'»„ 


Hb,(b,+bs)  -I... 
+ l^bf+BT)  *4 + 

Mr”-] 


_2A. 


f + 2B  f 
fl  + b2(b2+b4)*2  + 


“(^T  + b 


■ h20 


..EQ  (4) 

SEE  FIG. 2 

WHEN  C*0,V2f  CAN  BE  APPLIED  TO  A (AXISYMMETRIC) 
CYLINDRICAL  COORDINATE  SYSTEM, 

SAY  IN  R AND  Z : 


V f * A 


* 


0R 


£L 

0R 


.EQ  (5) 


FOR  A REGULAR  STAR  , THE  HARMONIC  OPERATOR 
BECOMES  (IN  A SIMILAR  MANNER  TO  EQ  (3)): 
hVfo  ■ [A  (f,  + »,)  +B  t f*  ♦ f<)  ♦ |Jl-  ( f,  - f4  ) 

- { A ♦ B ) 2 f0J  « h*D  ..  EQ  (6) 

SEE  FIG.  3 
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FOR  AN  IRREGULAR  STAR  ( ht*  M MM 
THE  HARMONIC  OPERATOR  BECOMES  (IN  A MANNER 
SIMILAR  TO  EQ  (4)): 


4 b3(b,*  b3)  + ttfb2+b4)f4* 

+ C b 1 b4 « _ ^ « 

R0  | b2(fc2  + b4)  2 b4(b2+b4)T4 


♦ 


_2_B _ C h / b2  ~b4  A a 

b2b4  R0  1 b2b4  T 


• h2  D ..  EQ  (7) 


SEE  FIG.  4 

EQS  (4)  AND  (7)  ARE  EMPLOYED  IN  THE 
PROGRAMMED  SOLUTIONS  FOR  CARTESIAN 
AND  CYLINDRICAL  COORDINATES,  RESPECTIVELY. 
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GRADED  MESH 


THE  MESH  SIZE  MAY  BE  REDUCED  IN  CRITICAL 
REGIONS , YIELDING  A HIGHER  RESOLUTION  WHERE  RE- 
QUIRED, WITHOUT  THE  COST  OF  AN  EXCESSIVE  NUMBER 
OF  NODES  OVER  THE  ENTIRE  DOMAIN  OF  THE  PROBLEM. 

THE  FINER  MESH  IS  TREATED  WITH  THE  SAME  EQUATIONS 
(EQS. (4)  & (7)  AS  THE  ORIGINAL,  BUT  WITH  THE  NEW 
SPACING,  h. 

BETWEEN  THE  ORIGINAL  (COARSE)  AND  NEW  (FINER) 
MESH,  HOWEVER,  THERE  EXISTS  AN  INTERMEDIATE  MESH 
OR  NET  THAT  REQUIRES  SPECIAL  TREATMENT.  THIS 
INTERMEDIATE  MESH  WILL  NOW  BE  CONSIDERED  FOR  BOTH 
CARTESIAN  (EQ  (4))  AND  CYLINDRICAL  (EQ  (7))  COOR- 
DINATE SYSTEMS.  NOTE  THAT  INTERMEDIATE  MESH  GRIDS 
ARE  "SQUARE".  THAT  IS,  ALL  ARMS  OF  THE  STAR  ARE 
EQUAL  (hi«h2-h3«h4-h) . 
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L v\;  + 

BECOMES 

(•^ ) [ f,  ♦ f,  4 f0j-  d/' 


DiL 

2 


.EQ (8) 


BECOMES 

.£i5>  + SlLu  # if  u / A±& Cl- 


( 


4R  ■»‘,-*V+(-t 


-<-4sfl-)4f  » D/Voh! 


WHERE 


' ‘9  Z 


4R 

• • EQ (9) 


><*♦«•> 


' 
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0 THINGS  TO  COME 
FLAT  PLATE  ANALYSIS 

The  deflection  (w)  of  a thin  plate  loaded  normal  to 
its  plane  is  described  by  the  fourth  order  partial  dif- 
ferential equation: 


a*w,  . d4W  84VJ 

3X4  ^9x2  8Y2  9Y4 


-g(x,y) 


Unfortunately,  because  of  the  many  plate  configura- 
tions possible,  a generalized  finite  difference  operator 
for  an  irregular  boundary  value  problem  perversely 
resists  formulation.  It  is  possible,  however,  to  replace 
the  fourth  order  equation  with  two  equations  of  the 
second  order,  since: 


V4f  x v2(V2f) 


The  two  equations  (see  eq.  12a, b)  are  of  the  same 
kind  as  that  obtained  for  a uniformly  stretched  and 
laterally  loaded  membrane  and  are  solved,  not  simul- 
taneously, but  sequentially.  The  input  or  forcing 
function  of  the  first  equation  is  the  lateral  loading 
g(x,y),  and  the  solution  variable  is  the  bending 
moment  vector  M at  each  node.  This  bending  moment 
vector  is  the  "loading"  input  to  the  second  equation 
which,  when  solved,  yields  the  deflection  of  the 
plate  at  each  node. 

This  solution  had  already  been  incorporated  into 
the  earlier  refresh  graphics  version  (CLYDE-274)  for 
simply  supported  plates  (M*0,w«0  on  the  boundaries), 
yielding  excellent  results.  Batch  and  storage  tube 
graphics  versions  are  planned,  tailored  to  solve  plate 
problems  with  a variety  of  edge  restraints  (simply 
supported,  built-in  or  "fixed",  free  edges,  etc..). 
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THE  BIHARMONIC  OPERATOR 


V4W- 


♦ 2 -2l)AL  ♦ '4 

9x4  dx^Y2  9Y4 


q 

D 


» 4-  -EQ  (10) 


CAN  BE  REPLACED  BY  TWO  SECOND  ORDER  EQUATIONS 

^2  » m s2tu  ■\2»u  . Q 


9x*  SY2  ax2  a Y2  * D 


EQ  (II) 


SINCE 

M* 

", 

Mx  + My  s 


*~o(*£L*MaY!^  * AND 
0( 

V BY2  ax2  ; 


-D(l  + M)( 


a2w  + a!w_ 


ax3 


aY! 


) 


INTRODUCING  A NEW  NOTATION 

l + p ° W SY2  ' 


EQ  (II)  MAY  BE  REPRESENTED  BY 

-q  • * EQ  ( 12a ) 


Q2M  + ?u_ 


ax2  by2 

a2w  a2w 
ax2  aY1 


EQ  (12  b) 


137 


TORSION  of  hollowed  shaft 


This  would  appear  to  be  a simple  matter  of  solving 
the  governing  PDE  over  a multiply-connected  boundary, 
were  it  not  for  the  uncertainty  concerning  boundary 
conditions.  The  actual  value  of  the  problem  variable  at 
the  boundary  was  not  important  in  the  torsion  application, 
only  the  difference  in  the  problem  variable  at  various 
points  mattered.  The  problem  variable  at  the  boundary 
could  be  assumed  to  have  any  value,  as  long  as  there  was 
only  one  boundary.  With  two  (or  more)  boundaries  the 
solution  calls  for  a different  approach. 

The  stress  function  is  obtained  as  the  superposition* 
of  two  solutions,  one  of  which  is  adjusted  by  a factor 
(k) . This  is  the  planned  programmed  solution  to  shafts 
with  a hole.  The  hole  may  be  of  any  shape,  size,  and 
location.  The  two  solutions,  to  be  combined,  are  shown 
at  the  right:  equations  and  boundary  conditions.  This 
capability  already  exists  in  CLYDE.  What  will  be  added 
is  the  solution  for  k and  then  the  final  superposition 
of  results.  Once  the  contour  integrals  are  taken  around 
the  inner  boundary  of  the  area  A_,  the  only  unknown, k, 
may  be  readily  obtained.  The  contour  integral  need  not 
be  evaluated  around  the  entire  boundary,  but  may  be 
taken  around  any  contour  that  encloses  that  boundary, 
and  includes  none  other  (for  example,  see  shaded  area 

V* 


i 

\ 


i 


1)  Reference  (5),  pgs  8,11,23 


I 
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THIS  P4CC  Dmtm  Mntmrmd) 


READ  INSTRUCTIONS  . I 


Q AUTHORS*  CONTROL  CARDS 

This  section  illustrates  additional  Control  Data 
SCOPE  control  cards  used  by  the  authors  to  catalog 
and  attach  their  input  and  output  data  files  and  to 
attach  the  executable  files  of  the  various  versions 
of  the  CLYDE  family.  This  is  for  illustrative  purposes 
only.  Users  are  enjoined  not  to  blindly  copy  these 
control  cards,  but  to  develop  their  own  file  names. 
Users  may,  and  in  fact  must,  exactly  duplicate  and  use 
the  ATTACH  coitmands  ending  in  "ID-MISDSEAD"  in  order 
to  use  the  CLYDE  programs.  They  should  not  reproduce 
the  authors'  COMMENT  cards  or  control  commands  con- 
taining "ID-ISAKOWER"  or  "ID-BARNAS".  Violators  will 
be  shot. 


AUTHOR'S  INITILIZATION  OF  TEKTRONIX  4014  TERMINAL  TO  RUN 
CLYDE-TEK  AT  360  BAUD. 


1.  Turn  terminal  (switch  1 foot  from  floor)  and  hard 
copier  on. 

2.  Put  LOCAL/LINE  switch  to  LOCAL. 

3.  Let  terminal  and  copier  warm  up. 

4.  Settings: 

CODE  EXPANDER  to  OFF 
CLEAR  WRITE  to  OFF 
ASB  toggle  switch  to  B 

ROTARY  BAUD  switch  to  300  (rear  of  terminal) 

ASCII/BCD  switch  to  ASCII  ( " " " ) 

Phone  SELECT  switch  CW  to  OPT 

5.  Hit  RESET  PAGE  key. 

6.  Hit  SHIFT  and  RESET  PAGE  keys  together. 

7.  Hit  SHIFT,  CTRL,  and  P keys  together. 

8.  Hit  RETURN  key  ( ©>. 

9.  Put  LOCAL/LINE  switch  to  LINE. 

10.  At  GRAY  phone,  depress  TALK  button. 

11.  Dial  361-6036  to  connect  to  CDC  6500. 

12.  When  computer  answers  with  constant  tone,  depress 
DATA  button  and  hang  up. 

13.  LOGIN  when  requested  by  system. 

Key  in  control  cards  as  required  and  run  program. 

All  user  keyboard  entries  are  followed  by  depressing  the 
RETURN  key  (this  is  symbolized  in  instructions  by  (r)). 

When  program  run  is  over  and  system  displays:  COMMAND- 

14.  Key  in  LOGOUT. 

15.  Hit  SHIFT  and  RESET  PAGE  keys  together. 

16.  Turn  off  equipment. 


JWf 


162 


300  BAUD  RUN  OF  CLYDE-TEK 


AUTHOR'S  INITIALIZATION  OF  TEKTRONIX  4014  TERMINAL  TO  RON 
CLYPE-TEK  at  4800  BAUD. 


1.  Turn  terminal  (switch  1 foot  from  floor)  and  hard 
copier  on. 

2.  Put  LOCAL/LINE  switch  to  LINE. 

3.  Let  terminal  and  copier  warm  up. 

4.  Settings: 

CODE  EXPANDER  to  ON 
CLEAR  WRITE  to  OFF 
AfcB  toggle  switch  to  A 

ROTARY  BAUD  switch  to  EXT  (rear  of  terminal) 
ASCII/BCD  switch  to  BCD  ( " " " ) 

PHONE  SELECT  switch  CCW  to  NORMAL 

5.  Hit  RESET  PAGE  key. 

6.  Hit  SHIFT  and  RESET  PAGE  keys  together. 

7.  At  GREEN  phone,  depress  TALK  button. 

8.  Dial  361-3785  to  convert  to  CDC  6500. 

9.  When  computer  answers  with  constant  tone,  depress 
DATA  button  and  hang  up. 

10.  LOGIN  when  requested  by  system. 

Key  in  control  cards  as  required  and  run  program. 
All  U8£r  entries  are  followed  by  depressing  the  RETURN 
key  (®). 

When  run  is  over  and  system  displays:  COMMAND  - 

11.  Key  in  LOGOUT. 

12.  Hit  SHIFT  and  RESET  PAGE  keys  together. 

13.  Turn  off  equipment. 
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(FILE  GENERATED  BY  CLYDE-BATCH) 
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CLYDE  ia  a computer  language  for  your  differential  equations. 
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29  (continued 


ere  are  three  versions  of  CLYDE.  This  document  describes 
the  capabilities  of  the  CDC  6000/TEKTRONIX  4014  storage  tube 
graphics  version  (CLYDE-TEK)  and  the  batch  version  ( CLYDE- B)  and 
also  serves  as  a preliminary  user's  manual.  An  earlier  version 
(CLYDE-274 ) , written  for  the  CDC  6500/1700/274  refresh  graphics 
facility,  is  described  in  MISD  Information  Report  73-1,  January 
1973,  and  includes  the  extension  of  the  solution  to  the  fourth 
order  stress  analysis  equation  for  flat  plates.  All  CLYDE 
versions  were  written  for  CDC  6000  host  computers  under  the  SCOPE 
operating  system  with  overlay  structures.  Two  applications 
covered  in  detail  in  this  document  are  steady  state  heat  conduct!) 
and  the  membrane  or  soap  film  analogy  of  torsion  of  bars  and 
shafts. 


